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We examine the QCD evolution of the helicity and transversity parton distribution functions 
when including also their dependence on transverse momentum. Using an appropriate definition of 
these polarized transverse momentum distributions (TMDs), we describe their dependence on the 
factorization scale and rapidity cutoff, which is essential for phenomenological applications. 
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I. INTRODUCTION 



Our understanding of the partonic structure of hadrons relies on the study of parton distribution functions (PDFs) 
and their extensions. In the last years, particular attention has been devoted to transverse-momentum-dependent 
parton distribution functions (TMDs). Standard coUinear PDFs are defined through collinear factorization theorems 
and obey the vifell-known DGLAP evolution equations [ll-Q- TMDs are defined through transverse-momentum- 
f~i \ dependent factorization and obey different evolution equations Here, for the first time we analyze these TMD 

Oh- evolution equations for two important distributions: the helicity and transversity TMDs. 

Factorization theorems are cornerstones of our understanding of hadron structure. They describe experimentally 
^ measured cross-sections in terms of perturbatively calculable hard parts and universal structures related to nonper- 
turbative parton dynamics, e.g., PDFs or TMDs. Factorization leads to well-defined evolution equations for the 
nonperturbative functions, which allows us to relate experimental measurements at different hard scales and perform 
' global analyses of PDFs, TMDs, and the corresponding fragmentation functions. 
J> ' The foundations of TMD factorization and evolution date back to Refs. B,H[ - However, important details related to 
0^ gauge invariance have been clarified only in the last decade (see, e.g., [1, (MtHI)- The first poof of TMD factorization 
was provided by Ji, Ma, and Yuan in Refs. [l^] while a complete definition of TMDs and rigorous proof of 
^""^ I factorization has been recently presented by Collins in Ref. @ and applied in Refs. [l3l - [l6j . Another definition has 
' been proposed in the context of Soft-Collinear Effective Theory by Echevarria, Idilbi, and Scimemi in Refs. [3 
CO . (see also the closely-related work of Cherednikov and Stefanis Refs. [ll, H^). The relation between the approach 
■ of Collins and that of Echevarria, Idilbi, and Scimemi has been analyzed in Ref. [2l| , with the conclusion that they 
' are essentially equivalent. In this work, we use the definition of Collins ^ and work in the framework of his TMD 
" . [ • factorization approach. 

^ ' The nonperturbative objects introduced in factorization theorems typically depend on a renormalization scale in 
the case of collinear PDFs, and a so-called rapidity cutoff C, in the case of TMDs. Physical quantities do not depend on 
these artificial scales, but only on the experimentally measurable hard scale of the process (e.g., the photon invariant 
, mass in a Drell-Yan process). 

For the description of a spin-1/2 target, eight independent TMDs can be introduced (at leading twist) [23 - |25| . At 
present, there exist explicit formulas for the evolution of two of them: the unpolarized TMD, fi{x, kx), and the Sivers 
function, f^{x, Ut). Here we consider the helicity distribution, gi{x, kr) and the transversity distribution hi{x, fe-r)- 
They are closely related to the collinear PDFs gi{x) and hi{x), whose collinear evolution is well known (see, e.g., 
Refs. Hiil). 

TMD evolution is related to transverse-momentum resummation and so called Collins-Soper-Sterman (CSS) for- 
malism [4]. For the polarized case of interest here, studies are presented in Refs. We will clarify in which sense 
our results correspond to the ones presented in the resummation literature. 
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II. TMD EVOLUTION 



We choose a particular process, namely Semi Inclusive Deep Inelastic Scattering (SIDIS). We denote with P and 
S the momentum and spin vector of the hadron target, and with the momentum of the detected hadron. With 
a single exchanged photon of momentum q, independent kinematic variables are: Q = a/— g^, x = Q^/2P ■ q, 
z = P ■ Ph/P ■ q, and the virtual photon's transverse momentum qt (in a hadron frame where the measured hadrons 
have zero transverse momentum). 

Details about TMD factorization and definitions are given in Refs. Here we summarize only the most 

important points. A SIDIS structure function in the form derived by Collins Q reads: 



Fuu,Tix, ql, g') - ^uu,t{Q-^ a*) j d^kr Spr /f (x, fer; A^, Cf) Dl{z, zpr; ^i, Co) 6^^'> (kr + qt - Pt) 



Yuu.T{Q,qT)+0{K/Q). 



(1) 



Here, 'Huu,t is the hard scattering part, fi{x,kT) is the TMD PDF for an unpolarized quark of flavor a in an 
unpolarised proton, and D1(z^zpt) is the unpolarized fragmentation function. The formula is similar to the parton- 
model expression (see, e.g., [25|), except for the dependence of the functions on fj, and the presence of higher-order 
terms in the hard scattering, and the presence of the correction term Y , which serves the purpose of correcting the 
expression of the structure function at large qx ~ Q, where the expression in terms of TMDs is not applicable. 
Analogous formulas can be derived for other structure functions containing other TMD PDFs and TMD FFs, see 
Refs. H^HOl- A denotes a generic hadronic scale, e.g., Aqcd or M. Note that we work in an approximation in 
which light cone momentum fractions that enter in the definition of Eq. ([1]) and usual Bjorken variable x are equal to 
each other. 

To correctly define the TMD PDFs and FFs in Eq. (HI), it is more convenient to define them in in transverse 
coordinate space (br-space) and then Fourier-transform the final result. The structure functions, as that in Eq. ([TJ, 
can be written as Fourier transforms of br-space expressions (see Sec. 2.2 of Ref. [soj). 

The proper definition of TMD PDFs requires the introduction of the so-called unsubtracted TMD PDFs together 
with further unsubtracted functions (sometimes called "soft factors" ) . Both of them contain rapidity divergences that 
are eventually canceled in the final definition of the the TMD PDFs. 

The unsubtracted TMD PDFs are defined as (dropping the flavor index) 



/— (x,6t;m,2/p 



Tr 



Vb) 
27r 



(P,5|V;(^/2)iy(C/2, 00,715(2/5)) 



t7 



W( 



e/2,oo,nB(yB))V'(-e/2)|P,5), (2) 



where = (0+,^^,6t), and we denote the functions with a tilde to indicate that they are defined in transverse 
coordinate space. The rapidity of the parent hadron is denoted by yp. An additional parameter yp is needed to 
regulate the light-cone divergences that result from using exactly light-like Wilson lines. This parameter is ultimately 
set to —00 in the final definition and the rapidity divergence is canceled by corresponding rapidity divergences in 
the soft function that we shall define below, see Ref. @. The subscript c indicates that only connected diagrams 
are included. The W{a, b; n) functions represent Wilson lines (gauge links) from a to 6 along the direction of the 
four- vector n. This direction is determined by the choice of the process It is essential for calculations of 

evolution for T-odd functions (Sivers and Boer-Mulders functions), but does not affect the discussion of the present 
paper, which focuses on T-even distributions. Light cone variables are defined such as — (a*' ± a^)/\/2 so that 
a ■ b = a^b^ + a~b'^ — ax ■ br- To obtain the expressions for the helicity and transversity TMD PDFs, the Dirac 
structure 7+ should be replaced by J~^j5, and 7~''7*75, respectively. 

For the proper definition of TMD PDFs, we need also to introduce an unsubtracted soft function that corresponds 
to the expectation value of a Wilson loop: 



S{o){bT\yA,yB) 



1 



(0|VF(bT/2, 00; nsY W{bTl2, 00; nA)W{~bTl2, 00; ns)VF(-6T/2, 00; n^)t |0) 



(3) 



In both ([2]) and ([3]), also transverse gauge links at infinity should be included [8]. However, when Feynman gauge 
is used their effects cancel in the final TMD PDF. Therefore we have not indicated the extra gauge links explicitly. 
The soft factor contains also self-interaction divergences that cancel in the final definition of TMDs. 

The complete definition of the TMD PDF in fey-space, given in Refs. @, is 



/i(x,bT;M,CF) = /r^"''(a;,6T;/^;yp- (-00)), 



S [o){bT; +00, Vs) 



\ S(Q){bT;+oo,-oo)S(^o){bT;ys,-oo) 



Zp Z2. 



(4) 
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Here, the "oo" arguments for the rapidity variables in the unsubtracted PDF and the soft factors are meant in the 
sense of a hmit. Ah field operators are unrenormalized, and Zp and Z2 are the PDF and field strength renormalization 
factors respectively. The soft factors on the right-hand side of Eq. (jl]) contain rapidity arguments ys ■ It is an arbitrary 
parameter which can be thought of as separating the plus and minus directions. 

To introduce polarization effects, we denote the Bloch 3-vcctor for a spin-1/2 particle moving in -\-z direction as 

p^{pT,\), (5) 

where A is the helicity and pr is transverse spin, so that in massless limit one has 



2 

If we define the following 4- vector 
then we can rewrite Eq. (|6]) as 



\^u{p)u{p) = \fil~ Xl5- X! TsPtTM • (6) 



P^ = (0+,0-,pt), (7) 



s 

this equation will prove useful when calculating Feynman diagrams. We also choose an appropriate Sudakov decom- 
position for vectors introducing two light cone vectors: 

p= (p+,0-,Ot) , (9) 
n= (0+,1-,Ot) . (10) 

The Fourier transform for TMD PDFs in D = 4 - 2e reads 

fix, br) = / d^-^'kTe-'''^"^f{x, kr) . (11) 



The inverse Fourier transform reads 

/(x, kr) ^ J |-^e'''-'^-/>, br) . (12) 

In order to address the problem of TMD evolution, we have first of all to compute the TMD PDFs in a parton-target 
model (see Fig. [TJa). We will study distribution of a parton of type j and momentum fc in a parton of type i and 
momentum p. At tree level (Fig. [T]a and the analogous case for gluons), the results for the TMD PDFs we want to 
consider are 

/— '"'''(^.fe^) - Sil-x)S'^'\kT)Si, (13) 
ffr"'""'(^,fcT,A) = S{l-x)5('HkT)Si, (14) 
/ir"""'-''(^,fcT,PT) = S{l~x)S(^HkT)5i , (15) 

where hi denote unpolarised distribution, helicity distribution and transversity distribution respectively. 

Including virtual gluon emission diagrams (see Fig. [1]) leads to multiplicative corrections to Eqs. ([T3l [T4l [TS)) of the 
form 

/r™''['"(a;, kr) = /^"''^(a;, fey) • C , (16) 
where C is the appropriate result for the virtual gluon loop: 

r ■ 2 2er ,j Tr(^(^-/)yi^) ^ .+u(2)ri. A n7A 
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(a) 




FIG. 1: Tree level diagram (a) and example of virtual gluon emission diagram (b) for the calculation of TMD PDFs in quark 
target model, p is the momentum of the parent quark containing quark k. 



The evolution of TMDs follows from their definitions, Eqs. Q. The rapidity evolution (with respect to C,f) is given 
by the Collins-Soper (CS) equation Q: 



ain/(a:,6T;M,CF)_^(^^.^^ (18) 



ainVcF 



where the function K{bT] fi) is defined as. 



K{bT■,^.)^l^\n (19) 
^ 2dys \S{bT;+oo,y,)J ^ ' 

Note that the rapidity evolution depends only on the Soft factor, which is independent of the polarization of the 
quark ^Ij. Note that the equation contains 5(6t) rather than S'(o)(bT)- Thus it is important to account for the UV 
renormalization factors ZpZ2 in Eq. (|4|. 

The dependence on the scale /i arises from renormalization group equations for both /(x, 6t; Mi C^") ^'^^ K(bT', /i). 
They are 

-^-lK\a{p)) (20) 



and 



din /i 

dlnf{x,bT;fJ,XF) 
dln/x 



lF(9{^l)■,CF/^^')■ (21) 



Here, g{ii) simply denotes the strong coupling with its dependence on the scale. The functions 'yK{g{fJ-)) and 
7f(5(m); Cp/pt^) are the anomalous dimensions of K{bT] /i) and f{x, br; /i, Cf) respectively. The fact that and 'yp 
are independent of bj is due to the fact that UV divergence arises from virtual gluon emission diagrams only, see 
discussion in Ref. j32{. As mentioned earlier, those diagrams give multiplicative factors, thus the overall result on 
evolution does not depend on the polarization of the quark. In conclusion, the results derived in Ref. [isj do not 
depend on the gamma matrix structure F used to define the specific polarized TMD, therefore they are universal for 
all polarization states and allow us to write down immediately the results for evolution of helicity and transversity 
TMD PDFs. 

There is a part where polarization is important. In the regime where kx is large compared to the hadronic 
scale, but still small compared to the hard scale (i.e., A <^ \kT\ <^ Q), TMDs can be calculated within a collinear 
factorization formalism |33l - l35| . This means that when bx is small but still larger than the inverse of the hard scale, 
i.e., 1/Q <C {brl ^ 1/A, Eq. (j?]) can be written as the convultion of a perturbatively calculable hard scattering 
coefficient and an integrated PDF: 

fiix, bT;fi, Cf) = ^ / —C,/,, {x/x, br; ^i, CF)fi (i; A*) + 0(A6t) , (22) 
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J J 
(a) (b) 

FIG. 2: Real giuon emission diagrams (a) and (b) of TMD PDF in quark target model, p is the momentum of the parent 
quark containing quark fc. Hermitean conjugate diagrams have also to be taken into account. 



the sum / goes over all quark and antiquark q, antiquark g flavors and gluon g. The functions fi{x;fi) are the 
ordinary integrated PDFs and the Cj/ji{x/x, bT] /i, Cf) are the hard coefficient functions. Similar expressions can be 
written for the helicity and transversity distribution. The hard coefficients will be different and will be denoted by 
AC* for helicity and dC for transversity. We will explicitly calculate them for helicity and transversity distribution 
function. 

Being independent of the type of initial hadron, the computation of the hard coefficients can be performed for the 
parton-target case [l^ . We can write perturbative results for the TMD distribution at small 67- up the the first 
order of perturbation expansion as (removing for convenience the dependence on the scales) 

/ W^' (x, br) = Cfl (x/x, br) ® ff}'' (x) + cf^^, (x/x, br) ® /H'^' {x) . (23) 

The symbol (E) means the convolution from Eq. ([2^ . The lowest order result for the hard coefficient is simply 

Cf/^, {x/x, 6t) - <5(1 - x/x)6j,, . (24) 

Using the results for the lowest order of /[''I from Egs. p^ [Til [T5|) integrated over kx we obtain 

cfl^ix,bT) = f7J\x,bT)-f['^'ix). (25) 

This expression represents the recipe to compute the hard coefficients at order as- Analogous formulas hold for the 
hard coefficients AC and 5C of the helicity and transversity distributions. 

For the TMD PDF of a quark in a quark we follow the steps of Refs. [1, [l^. To deal with eikonal propagators, we 
use the Feynman rules from Refs. [6, 33]. From diagrams (a) and (b) in Fig. [2] we obtain 



f[]^\x,kr)l^g'fi'^Cp6i I 



dk-d'^-^H Tr {^ip - Dlcf (1 - A75 + 75^^,) r){t - /)) 
(27r)3-2^ 4((^~~02T^i)2 (26) 

X 6^^\It + kT)5[p+ -1+ - k+)5{e) , 

dk-d'^-^H {up - t)if (1 - + 75^j.)) 

(27r)3-2<= 4((p-/)2 + ze)(n-/-ie) (27) 
X (5(2)(Zj. + fcT)'5(p+ - /+ - fc+)(5(P) . 



The Dirac structure iji has to be replaced by ^^75 and ^^7*75 for the helicity and transversity distributions, respectively. 

Calculations of integrals in Eqs. (I26I27|) go along the lines of Ref. [6]. After estimating the trace in 4 — 2e dimensions 
we evaluate the integral over dl~ first by closing the contour at infinity and using Gauchy's integral theorem. Then 
we evaluate the reminding dV^ and c?^~^'^Zt integrals utilizing delta functions. Finally we compute MS counterterms 
by prescription from Ref. 

Using Eq. (P5|) and results of Eqs. (|26I27|) along with the soft subtraction factors (see Appendix A of Ref. [l^) for 
the TMD PDF for finding a quark of flavor ]' in a quark of flavor j we flnd to order a^. 
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S{1 - x) 



In' 



1 



X 



1 -X 



fibj 



In 



2(1--)+ 



In 



+ Oial), (28) 



In 



5{1 - x) 



/2e- 



7E 



In" 



In 



+ 2(1-^) + 



^62 



+ (29) 



SCj'/j{x,bT; h-Xf/iJ' ) =(5j'j5(l-x) +(5 



2.T 

1 - x 



5(1 - x) 



-In' 



fibr 



In 



2e-n^ /C. 



^67 



+ 0(«s)- (30) 



for unpolarised, helicity, and transversity TMDs respectively. The strong coupling 0:3 is evaluated at a scale fi, and 
the number of active flavors is e Nf. The usual SU{Nc) color factors are Cp = (N^ - l)/(2A^c), ?> = 1/2. 





FIG. 3: Quark in a gluon diagram for TMD PDF, p is the momentum of the parent gluon containing quark k. 

The calculation for quark in gluon follow the lines of Refs. 0, [l3|. We utilize the diagram from Fig. [3l Note that 
at the order there is no contribution from soft factor subtraction and all the diagrams with gluon attached to the 
Wilson line give zero if an appropriate choice of gluon polarization is made, see Refs. @. Furthermore, the diagram 
from Fig. [3] gives zero for transversity, since using the transversity projector 7^7*75 the Dirac trace contains odd 
number of Dirac matrices. This result is well known in the literature. 

Using Fcynman rules from Ref. ,6:] we obtain for diagrams in Fig. [3] 



T ^2, ,2c 



4(1 -ej 
4(1 -3e)' 



(27r) 



4-2e 



(fc2 



- le) 



s{{p-ky 



dk-d^-^H Tr {iToht.{$ - fh.^) 



(fc2 + ieY 



5{{p-kf 



(31) 
(32) 



for unpolarised and helicity distributions accordingly. Here 



1 



nV) 



(33) 
(34) 
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In conclusion, for the gluon contributions we obtain 

C,/, (x, br; /i, (p/f^') = ^ |ln (^) + (1 - ^)') + - x)| + 0(a2). (35) 

AC,/, (a;, 6t; A^, C^^/m') = ^ jln (^) (2^^ - 1) + (1 - ^)| + ©(a^- (36) 

SCg^jix, br; 11, Cf/a^') = + 0(a2). (37) 

for unpolarised, helicity, and transversity TMDs, respectively. Note that there are no contributions from Soft factor 
at this order to quark in a gluon coefficient functions. 

Eqs. (1^150)) and Eqs. ([Ml [57]) represent the original results of this paper. They allow us to write epxressions for the 
helicity and transversity TMDs that fulfill TMD evolution equations, Eqs. ([T51 120[ [^T|) . and have a behavior at high 
transverse momentum that matches perturbative calculations. The solution for a TMD for flavor i can be written in 
a compact way as (see Refs. [1, HI]) 

Jl [x, bT-^ti, Cf) = Y.(^f/^ ® fl) M6)e^^''*^^''''''^"^e'"^''"^'" ^/Sp(a;, br) (38) 



Analogous formulas hold for the helicity, gi, or transversity, hi, distributions. The sum goes over all quark and 
antiquark flavors and include also gluon, j — q, q, g. Appropriate coefficient functions should be used in each case from 
Eqs. (PSI (E) denotes the convolution in longitudinal momentum fractions of Eq. (1^^ . The scale fib is chosen 

appropriately to ensure the optimal convergence of perturbative series. The function /np denotes the nonperturbative 
part of the TMD and has to be fitted to experimental data. In the literature, it is usually parametrized as a Gaussian, 
although there is no fundamental reason for this choice. 

In order to be able to use Eq. (|22p also at large br, the so-called 6* prescription can be introduced. The function 
6* serves the purpose of freezing the value of bx, preventing it from becoming larger than a certain value and avoid 
regions where the perturbative running coupling as becomes divergent (at very large values of bx, or equivalently, at 
very small transverse momentum). A common choice is to set 



6* = 



,2 ' 
max 



(39) 



but other functional forms can be explored, as well as different prescriptions (e.g., the complex-6 prescription of 
Ref. [l^d^l)- Any change in the prescription should also be combined with a change of the nonperturbative function 

/i 



NP- 



The perturbatively calculable function S'(6*; /x^, /x, Ci?) reads 



S{b^;fj.b,fJ;Ci- 



In 



fib 



lF{g{^l')■,l) 



In- 



-T-7k(5(m')) 



(40) 



and expressions for K,^f, and at order as can be found in Appendix B of Ref. [Li 
We adopt also these other choices ^ 



fib = 2e = feo/^*, 



Cf = Q^ 



Cfo — Qo- 



(41) 



Different options can be explored in order to test the sensitivity of the final results to the scale choice. Note that 
K{b^;fib) = at this order with this choice. 



Using Eq. (|4T|) we find to order as for Eqs. (|28l29l30p 

Cj, /j{x,b*;fib) ^5j'j5{l~ 

ACj,/j{x,b^,;fib) = Sj'j6{l - 
SCj,/j{x,b^;fib) = (5j'j(5(l - x)+0{ai) 



x) + 5,.^ 



27r 
asGp_ 
27r 



(1 



0{al), 
0{a's): 



(42) 

(43) 
(44) 



^ Note that parameter bo is dimensionless, so that feo/ii* has the dimensions of energy, GeV. 
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and for Eas. (l35l36l37l) . 



^:E(l-x)+0(a2) 

TT 
TT 



for unpolarised, helicity, and transversity TMDs, respectively. One can see that C. 



j'/j 



(45) 

(46) 
(47) 

ACj/ /j and the difference 



between Cji /j and SC^i is the absence of a., contribution. 



At this point, a discussion on the comparision with the CSS hterature is in order. The restructuring of the formahsm 
in terms of TMD definitions makes it non-trivial to map them onto the classic CSS components. However, we can 
check that the final result for the structure functions match. For unpolarized DIS, results in the CSS formalism are 
presented, e.g., in Ref. [s^. Polarized DIS has been discussed in Ref. For instance, we can compare the results 
for the unpolarized structure function of Eq. ([1]). We need the expression of the hard scattering, which has been 
reported in Ref. [1^ 



1 



rin -;in^ -4 



(48) 



If we insert the expression for the unpolarized TMD, Eq. the analogous expression for the fragmentation 

function, Eq. (31) of Ref. and the hard scattering, Eq. (^5]) . in the formula of the structure fuction, Eq. ^ 
and we adopt the choices of Eg. (PTI) . we recover the standard CSS results for DIS at order as (see, e.g., Eq. (45) of 
Ref. [as'l, or Eq. (32) of Ref. [29^). Other structure functions for polarized DIS require the study of different TMD 
PDFs and FFs, together with a careful assesment of the possibility to match the low and high transverse momentum 
results (see, e.g., Ref. [35l.l40l. [4l]|'). Since here we are concerned only with helicity and transversity TMDs, results can 
be compared with those of Ref. f2!l| . 

Our results can be compared also to the Drell-Yan case. The unpolarized case has been discussed in many papers, 
e.g., Refs. i,|42-lil)- In this case, however, the hard scattering of Eq. (^51 has an extra 7r^/2 in the square brackets, 
as in Eq. (6) of Ref. [39|). Our results can be compared to those for doubly longitudinally and transversely polarized 
Drell-Yann [ii,!!^. 

Apart from the full structure function, we can focus our attention on the coefficient functions: our results expressed 
in Eqs. dH US IMl) and in Eqs. (0211111171) correspond to Eq. (37, 38, 41, 42, 43) of Ref . [29|]. Differences occur only 
for the terms with (5(1 — x) and can be ascribed to the contribution from the hard scattering. These contributions turn 
out to be the same for all three cases, which indicates that the hard scattering is the same for the relevant structure 
functions. This has been already observed, albeit with a different definition of the hard scattering, in Refs. (12i. i4l|. 
Since the hard scattering is different in different processes, in the CSS formalism the coefficients are different in DIS 
and Drell-Yan scattering. This can be checked by comparing, e.g., Eq. (38, 42) of Ref. |29[ with the coefficient Acg 
in Eq. (15) of Ref. [46] and Eq. (38, 43) of Ref. ^ with the as part of Eq. (8) of Ref. The difference between 
these results is due to the factor 7r^/2 that we discussed previously. 



III. PHENOMENOLOGY 



Let us finally present results of TMD evolution using a test function at initial scale Qo = 3.2 GeV. We choose the 
coUinear functions in the following form: 

a;uo(a;) =a;^^o(a;) = .x"■^(l-2:)"■^ xuoix) = xdo{x) = 0, xgo{x) = x°-^ {I - x)"-^ . (49) 

Note that the choice is arbitrary as our goal is just to demonstrate the results of evolution. These initial distributions 
will be assumed the same for unpolarised (/i), helicity (51), and transversity (hi) distributions. Results with realistic 
initial functions will be presented elsewhere. 

Evolution in br is identical for all three functions, see Eq. (|38|. The convolution in Eq. ((38l) takes three different 
forms. For the unpolarised distribution, using Eqs. (I42[ I45p we find 

{Cj/^<E>fi){x,b,;^ib) = fi{x,fib) + ^ / —(l-^)fi{x,fib) + ^^ / —(l-^)g{x,fib), (50) 
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X kj^ (GeV) 
fi ^ 

FIG. 4: Evolution of P {x, kT;Q, Q^) and P{x;Q,Q'^) at two different scales Q = 3.2 GeV (two upper plots), fO GeV (two 
middle plots), and fOO GeV (two bottom plots). Solid line corresponds to unpolarised evolution, dashed line corresponds to 
helicity evolution, and dotted line corresponds to transversity evolution. 



where g{x, ^b) corresponds to the gluon distribution. 

For the hehcity distribution, using Eqs. we obtain 

{ACj/,(E)gi)[x,K;fib) = gi{x,fib) + ^ — / — [1 - - j gi{x, fib) + / — [1 - -j Ag{x, fib) ■ (51) 

Finally, for the transversity distribution, using Eqs. (|44[ 1471) we obtain 

^{SCj/, «) h{){x,b^;^ib) = hi{x,fib) ■ (52) 

3 

Note that, as is well known for CSS resummation, at this order TMD evolution does not have mixing of different 
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flavors and no mixing with gluon TMD, even though colhnear PDFs mix as can be seen from Eqs. (j50l51l52p . This 
means also that the gluon TMD cannot be studied via scaling violations at least at this order. 

We perform DGLAP evolution for the collinear functions fi{x, iib),gi{x, fib), hi{x, /i^) using the HOPPET evolution 
package [isj . 

The choices of non-perturbative functions that enter in Eq. ((38|) are the following (we use again the same function 
for all three polarization cases just for illustration purposes): 

/^p(x,6t) = exp (^JlBl^ , g^^b^) = -gh%, (53) 

where {k\) ^ 0.25 (GeV^), g = 0.2 (GeV^). We also choose 6™ax = 1-5 (GeV"!). 

The choice of /^p corresponds to non-perturbative functions used in analysis of TMD functions at tree-level by 
the Torino-Cagliari-JLab group [i^ [s^. The choice of gKibr) is motivated by the so-called LBNY fit of Drell-Yan 
cross-sections using CSS resummation formalism ^51]. 

We will show the evolution of the TMD functions, as well as their integral over fc^ up to the value of Q, which we 
conventionally referred to as their 0*'' fc^-moment: 

f'^{x,kT; hXf) , Ci^) = 27r / kTdkTfq/pix,kT; h,Cf) ■ (54) 

Jo 

Note that this 0**^ fcT-moment af a TMD functions /^(a;; Cf) should not be confused with collinear PDF /^(a;; /i). 

In Fig. (U]) we show results of the evolution of /^(a;, fcr; Q, Q^) and /"^{x; Q, Q^) at three different scales Q — 3.2, 10 
and 100 GeV, for the unpolarized, helicity and transversity distributions. As one can see from Fig. after evolution 
the three functions become wider and are still very similar. The only appreciable differences are at higher transvese 
momentum. The O"' /cT-moments after evolution show some differences at low x. Compare our results with results of 
collinear evolution of gi and hi presented, for example in Ref. [52l|. As in the case of collinear evolution, also in TMD 
evolution hi becomes smaller than gi under evolution and the difference grows with Q. The reason is the absence of 
al contributions to coefficient functions of transversity. 

We also note that TMD functions at initial scale are very much similar to TMD functions parametrized at tree level 
[49I [50| which thus justify extraction of those functions at tree level. Observables at different characteristic scales 
however should be described using TMD evolution. 

One can also observe that the so called Soffer bound [53l|: \hi{x, Q^)\ < \ {fi{x, Q^) + gi{x, Q^)) , is also satisfied 
for TMD distributions f'^{x,kT;Q,Q^) and f''{x;Q,Q^) numerically. Let us remind that Soffer bound for collinear 
densities was shown to be preserved at LO accuracy in Ref. [1^ and at NLO accuracy in Ref. [23] ■ We set aside the 
discussion of Soffer bound for TMD functions for a separate publication. 

IV. CONCLUSIONS 

In this paper we calculated the evolution of the transverse-momentum-dependent (TMD) helicity and transversity 
distribution functions. We adopted the definition of TMD PDFs as given by Collins in Ref. [6,]. We provided explicit 
formulas for all coefficient functions at as- The results of this paper can be readily used in TMD phenomenology. 

As an illustration, we calculated the unpolarized, helicity and transversity TMD distributions at different scales, 
starting from the same initial conditions. The final results are very similar. Their 0'^ fc^-moments differ at low x. We 
observed that if started from equal initial conditions, helicity TMD distribution gi becomes smaller than unpolarised 
/i distribution and transversity hi becomes smaller than helicity gi TMD. 
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